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1. INTRODUCTION 

The idea that our Universe be a 3 + 1 dimensional membrane in a higher dimensional 
'bulk' spacetime has opened new exciting prospects for cosmology, for reviews see Oi 
0]. In the simplest braneworlds motivated by string theory, the standard model particles 
are confined to the brane and only the graviton can propagate in the bulk. Of particular 
interest is the Randall- Sundrum (RS) model Jsl, 13], where the bulk is 5-dimensional 
anti-de Sitter space, AdSs. If the so called RS fine tuning condition is satisfied, it can be 
shown that gravity on the brane 'looks 4-dimensionar at low energies. 

Within this model, cosmological evolution can be interpreted as the motion of the 
physical brane, i.e. our Universe, through the 5d bulk. Such a time-dependent boundary 
does in general lead to particle production via the dynamical Casimir effect 

Of course one can always choose coordinates with respect to which the brane is 
at rest, e.g. Gaussian normal coordinates. But then usually (except in the case of de 
Sitter expansion on the brane [6J), the perturbation equation describing the evolution of 
gravitons is not separable and can be treated only with numerical simulations 10, [sl,^]. 
Furthermore, in a time-dependent bulk a mode decomposition is in general ambiguous 
and one cannot split the field in a zero mode and Kaluza- Klein (KK) modes in a unique 
way. 

Based on the picture of a moving brane in AdSs, we have studied graviton production 
in an ekpyrotic type scenario ifioll where our Universe first approaches a second static 
brane. After a 'collision' the physical brane reverses direction and moves away from 
the static brane, see Fig. [TJ For an observer on the brane, the first phase corresponds 
to a contracting Universe and the collision represents the 'Big Bang' after which the 
Universe starts expanding. 



Here I report on the results which we have obtained in our previous papers [lllLll2l . ll3|l . 



We have found that the energy density of KK gravitons in AdSs scales like stiff matter. 
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FIGURE 1. Two branes in an AdSs spacetime. The physical brane is to the left. While it is approaching 
the static brane its scale factor is decreasing, the Universe is contracting, and when it moves away from 
the static brane the Universe is expanding. The AdS curvature radius L (dashed line) and value of the scale 
factor of the brane metric as function of the extra dimension y (light (blue) line) are also indicated. 



Pkk °^ ci^^, here a denotes the scale factor defined in Eq. Therefore, KK gravitons 
in AdSs cannot represent the dark matter in the Universe. This finding is in contrast with 
the results of Ref. [14] and we comment on this below. We have also found that in the 
early Universe the back reaction from KK gravitons on the bulk geometry is likely to be 
important. 

Finally, we have derived a limit for the maximal brane velocity, the bounce velocity, 
Vh^O.l in order not to over-produce zero-mode (i.e. 4d) gravitons, the energy density 
of which is constrained by the nucleosynthesis bound. We have calculated the spectra of 
both, the zero-mode and the KK gravitons. In Refs. [[I1IIT2I] we have, however, neglected 
a term linear in the brane velocity v in the boundary conditions. In our latest work, 
Ref. [|l3|] we derived a method which includes this term and allows to treat the problem 
without any low velocity approximation. We have shown that the low velocity results 
previously obtained are not modified. 

The remainder of this paper is organized as follows. In the next section we present 
the basic equations for the evolution of tensor perturbations (gravitons) and we explain 
why it is not straight forward to include the velocity term of the boundary condition. In 
Section 3 we quantize the system. In Section 4 we discuss our results and in Section 5 
we conclude. 



2. A MOVING BRANE IN AdSs 
2.1. The background 

In Poincare coordinates (x^) = (?,x,y) with x = (x^.x^.x^) and A = 0, ...,4, the AdSs 
(bulk) metric is given by 

ds^ = gABdx^dx^ = — [-dt^ + dijdx'dx^ + dy'^] , (1) 

where /, j = 1,2,3 and L is the AdSs curvature radius which is related to the bulk 
cosmological constant by the 5d Einstein equation, —A = 6/L^. The physical brane 
representing our (spatially flat) Universe is located at some time dependent position 
3; = yb{t) in the bulk, and the metric induced on the brane is the Friedman-Robertson- 
Walker metric, 

ds'^ = a^{ri)[-dr]^ + dijdx'dx^] , (2) 
with scale factor aij]) which is given by the brane position, 

<r]) = . (3) 

The conformal time 7] of an observer on the brane, is related to the bulk time t via 

dr] = \/ 1 - v^dt = Y^^dt . (4) 
Here we have introduced the brane velocity 

dyb LH 1 

V = — = , and r = , . (5) 



H is the usual Hubble parameter 

1 da 
dr\ 



H=-^^ = a~^J^ = -L-^rv- (6) 



The brane dynamics, as a result of the second junction condition, is determined by the 
modified Friedmann equation yj 



where o is the brane tension, p the energy density on the brane, and we assume the RS 
fine tuning condition [3] 



We define the string and Planck scales by 
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(9) 



Note that the RS fine-tuning condition is equivalent to 

K-5 = K4L or — = 



(10) 



2.2. Tensor perturbations 

We now consider 3d tensor perturbations hij{t,x,y) of the spatial three-dimensional 
geometry on this background. The perturbed bulk metric reads 



ds^ = ^ [-dt^ + (dij + 2hij)dx'dx-' + dy^] . 



(11) 



Tensor modes satisfy the traceless and transverse conditions, h'- = dih'j = 0. These con- 
ditions imply that hij has only two independent degrees of freedom, the two polarization 
states • = X , -|-. We decompose hjj into spatial Fourier modes. 



hij{t,x,y) 



d^k 



e"'-e-,(k)/z.(r,3;;k) 
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(12) 



where e*^(k) are unitary constant transverse-traceless polarization tensors which form 
a basis of the two polarization states • = x,-f . Since we assume parity symmetry, we 
shall neglect in the following the distinction between the two graviton polarizations and 
consider only one of them. We then have to multiply the final results for e.g. particle 
number or energy density by a factor of two to account for both polarizations. 

The perturbed Einstein equations and the second junction condition lead to the fol- 
lowing boundary value problem 



y y 



/i(f,y;k) =0 inthebulk, ]^ 



and 



Y{vdt + dy)h 



yb{t) 



on the brane . 



(13) 



(14) 



We introduce also a second, static brane at position ys, which requires the additional 
boundary condition 



on the static brane 



(15) 



Eq. (1131) is the Klein-Gordon equation for a minimally coupled massless mode in 
AdSs , i.e. the operator acting on h is just the Klein-Gordon operator 



1 



(16) 



Equation (fT4l) is the time-dependent boundary condition (BC) coming from the fact 
that the moving brane acts like a "moving mirror" for the gravitational perturbations. 
Only in the rest-frame of the brane do we have pure Neumann BC. In a generic frame 
we have the Lorentz transformed BC which contains a velocity term vdf. 

We assume that the brane is filled with a perfect fluid such that there are no anisotropic 
stress perturbations in the brane energy momentum tensor, i.e. there is no coupling of 
gravitational waves to matter. If this were the case, the r.h.s. of Eq. (fT4l) would not be 
zero but a term coupling hij to the matter on the brane, see Eq. (2.25) of [12]. 

For the tensor perturbations the gravitational action up to second order in the pertur- 
bations reads 



2k5 J J Jy,{t) L 



(17) 



One factor of two in the action is due to Z2 symmetry while a second factor comes from 
the two polarizations. 



2.3. Dynamical Casimir effect approach 



The wave equation (1131) itself has no time dependence and simply describes the 
propagation of free modes. It is the time dependence of the BC (fT4l) that sources the 
non-trivial time-evolution of the perturbations. As it is well known, such a system of a 
wave equation and time-dependent BC lead, within a quantum mechanical formulation, 
to particle production from vacuum fluctuations. In the context of the photon field 
perturbed by a moving mirror this goes under the name "dynamical Casimir effect" jsj]. 

In 1I12I1 we have extended a formalism which has been successfully employed for the 



numerical investigation of photon production in dynamical cavities [|15L I16L 1171 ] to the 
RS braneworld scenario. We have studied graviton production by a moving brane, which 
we call dynamical Casimir effect for gravitons, for a bouncing braneworld scenario. 

However, in order to solve the problem, we have neglected the velocity term in 
Eq. (Ull). The ansatz 

h = Y,aa{t)e''"'"'Mt,y) , Col = k^ + ma{tf 

a 

then leads to a Sturm-Liouville problem for the instantaneous eigenfunctions 0c<(^j) 
which satisfy 

-d^ + -dy](^a=ml(^a ■ (18) 



y 



The solutions of (fTSi) are 



m = ^M^. (19) 

^y^-yli^) 

^n{t,y) = Nn{t)y^C2{mn{t),yb{t),y) with 
Cv(m,jc,j) = Yi{mx)Jv{my)—Ji{mx)Yv{my). (20) 



The function 0o is the zero mode which corresponds to the ordinary (3 + l)d graviton on 
the brane while the ^„ are the KK modes. The masses m„ are determined by the boundary 



condition at the static brane, see, e.g. [|18ll for more details. Since satisfies Neumann 
boundary conditions, we know that the solutions {^a)a form a complete orthonormal 
set of functions on the interval [yz7(?)5 Ji] normalized by the scalar product 



Jyh(t) y 



yh{t) 

Therefore, any general solution which satisfies Neumann BC can be expanded in these 
instantaneous eigenfunctions. If we add the term vdt to the boundary condition this 
feature is lost, and we can no longer expect to find a complete set of instantaneous 
eigenfunctions. 

However, since the entire effect disappears when the velocity tends to zero, neglecting 
a term which is first order in the velocity seems not to be consistent. This problem led us 
to search for another approach which is discussed in Ref. where we transform to a 
coordinate system where the velocity term disappears identically. There also show that 
for low velocities v < 0.3, say the corrections obtained with this consistent treatment are 
below a few percent. We therefore ignore it in the following. 



3. QUANTIZATION 
3.1. Equation of motion 

The gravitational wave amplitude h{t,y;'k) subject to Neumann boundary conditions 
can be expanded as 

h{t,y;k) = a/ yf £ qa,k{t)Mt,y) ■ (21) 

The coefficients <?a.k(0 are canonical variables describing the time evolution of the 
perturbations and the factor ^/k^JI? has been introduced in order to render the qa.k^ 
canonically normalized. For h{t,y,x) to be real, we have to impose the following reality 
condition on the canonical variables, 

= ^a-k- (22) 

One could now insert the expansion (|2T]) into the wave equation (fT3] ), multiply it 
by ^p{t,y) and integrate out the dependence by using the orthonormality to derive 
the equations of motion for the variables qa,k- However, as we explain in Refs. [1T2L 
[13I1 . a Neumann boundary condition at a moving brane is not compatible with a free 
wave equation. The only consistent way to implement Neumann boundary conditions is 
therefore to consider the action (flTl) of the perturbations as the starting point to derive 
the equations of motion for qaM- Inserting (|2T]) into ([17]) leads to the action 



E [^ccP (^a,k4j8.-k + ^a,-k4j6,k) + A^ajS^a.k?^ -k] \ ■ (23) 
aj8 ^ 



We have introduced the time-dependent frequency of a graviton mode 

03l.k = \Jk^ + nil, (24) 
and the time-dependent coupling matrices 

Mafi = (Ma^^p) , (25) 
Nap = {dt^a,Mp)=Y.^ayMpy={MM^)ap, (26) 

y 

which are given explicitely in Ref. fl2l] (see also fis']). The equations of motion for the 
canonical variables are the Euler-Lagrange equations from the action (|23l) . 

'?a.k + «a,^^a,k + £ [^^a -Map] 4/3 .k + £ - ^jS.k = • (27) 

The motion of the brane through the bulk, i.e. the expansion of the universe, is en- 
coded in the time-dependent coupling matrices M(^^ and Nap ■ These mode couplings 
are caused by the time-dependent boundary condition dyh»{t,y) = which forces the 
eigenfunctions 0a(?,j) to be explicitly time-dependent. In addition, the frequency of 
the KK modes 0)a.k is also time-dependent since the distance between the two branes 
changes when the brane is in motion. Both time dependencies can lead to the amplifica- 
tion of tensor perturbations and, within a quantum treatment which is developed below, 
to graviton production from vacuum. 

Because of translational invariance with respect to the directions parallel to the brane, 
modes with different k do not couple in (|271) . The three-momentum k enters the equation 
of motion for the perturbation only via the frequency cOa.k- Equation (ITTI) is similar to 
the equation describing the time evolution of electromagnetic field modes within a three- 
dimensional dynamical cavity [[l6|] and may effectively be described by a massive scalar 
field on a time-dependent interval (vfj. For the electromagnetic field, the dynamics of 
the cavity, or more precisely the motion of one of its walls, leads to photon creation 
from vacuum fluctuations. This phenomenon is usually referred to as dynamical Casimir 
effect. Inspired by this, we call the production of gravitons by the moving brane the 
dynamical Casimir effect for gravitons. 



3.2. Quantization 

Asymptotically, i.e. for t ±00, the physical brane approaches the Cauchy horizon 
(jb 0)5 moving very slowly. Then, the coupling matrices vanish and the KK masses 
become constant, 

lim M(y,B{t)=0 , lim ma(0= const. Va,/3 . (28) 



In this limit, the system (1271) reduces to an infinite set of uncoupled harmonic oscillators. 
This allows to introduce an unambiguous and meaningful particle concept, i.e. the notion 
of (massive) gravitons. 

Canonical quantization of the gravity wave amplitude is performed by replacing the 
canonical variables qa.k by the corresponding operators qa,k 



h{t,y;k) = J ^ J^<?a,k(0<^>a(^j) 



(29) 



Adopting the Heisenberg picture to describe the quantum time evolution, it follows that 
4a,k satisfies the same equation (l27l) as the canonical variable ^a,k- 

Under the assumptions outlined above, the operator 4a,k can be written for very early 
times, t < ?in, as 



1 



^in „-'(Oa^kt I -int JO)]^i,t 



+ a 



(30) 



where we have introduced the reference frequency 



(31) 



This expansion ensures that Eq. (I22l) is satisfied. The set of annihilation and creation 
operators {a'^^, k) corresponding to the notion of gravitons for t < is subject to 
the usual commutation relations 



^in -int 



-m -in 
"a,k5"a'.k' 



"a,k'"a',k' 



0. 



(32) 
(33) 



For very late times, t > fout, i-C- after the motion of the brane has ceased, the operator 
qa,k can be expanded in a similar manner. 



out; 



^"ce,-k^ 



'a,k<= 



with final state frequency 



< ^ = ^a,k{t > tout) 



(34) 



(35) 



The annihilation and creation operators {«a"k5^a"k^} correspond to a meaningful defi- 
nition of final state gravitons (they are associated with positive and negative frequency 
solutions for t > tout) and satisfy the same commutation relations as the initial state op- 
erator^. 



Of course the brane never really stops moving, but before a certain time fin and after a certain time fout 
the motion is so slow that no particle production takes place. We have chosen these times sufficiently early 
(rsp. late) so that the numerical results are independent of their choice. 



Initial |Ojin) = |0,f < fin) and final |0,out) = |0,f > fout) vacuum states are uniquely 
defined viaEI 

«a,k|0.in) = , a°"|,|0,out) = , V a, k . (36) 

The operators counting the number of particles defined with respect to the initial and 
final vacuum state, respectively, are 

ju'm _-int-in iCrout _ -outf -out /o^n 
^^a,k — '^a,k^a,k ' ''^a,k — ^a,k ^a,k • ' ) 

The number of gravitons created during the motion of the brane for each momentum k 
and quantum number a is given by the expectation value of the number operator iV°"^ 
of final-state gravitons with respect to the initial vacuum state |0, in): 

^at=MN^V^M- (38) 

If the brane undergoes a non-trivial dynamics between ?out we have «^"^|0, in) 7^ 

in general, i.e. graviton production from vacuum fluctuations takes place. 



4. RESULTS 

4.1. Energy density 

For a usual four-dimensional tensor perturbation /i^v on a background metric g^v an 
associated effective energy momentum tensor can be defined unambiguously by 

T^iv = ^{hap\\^h''^ny). (39) 

where the bracket stands for averaging over several periods of the wave and "||" denotes 
the covariant derivative with respect to the unperturbed background metric. The energy 
density of gravity waves is the 00-component of the effective energy momentum ten- 
sor. We shall use the same effective energy momentum tensor to calculate the energy 
density corresponding to the four-dimensional spin-2 graviton component of the five- 
dimensional tensor perturbation on the brane, i.e. for the perturbation hij{t,x,yb). For 
this it is important to remember that in our low energy approach, and in particular at 
very late times for which we want to calculate the energy density, the conformal time 
rj on the brane is identical to the conformal bulk time t. The energy density of four- 
dimensional spin-2 gravitons on the brane produced during the brane motion is then 
given by 

P = (^(o^n\hij{t,x,yb)h\t,x,yh)\0,m^Y (40) 

Here the outer bracket denotes averaging over several oscillations, which we embrace 
from the very beginning. The factor l/a^ comes from the fact that an over-dot indicates 



Note that the notations |0,f < fin) and |0,f > fout) do not mean that the states are time-dependent; states 
do not evolve in the Heisenberg picture. 



the derivative with respect to conformal time t r]. The detailed calculation given in 
Ref. iilfl leads to 

2 f d^k 

P = ^I^J J^^a,k^oc,km^{a) (41) 

where again yya,k{t) is the instantaneous particle number and % is related to value of 
the wave function on the brane by 

The factor two reflects the two polarizations. At late times, t > ?out, after particle creation 
has ceased, the energy density is 

P = ^ I / «S <k (a) . (42) 

This expression looks at first sight very similar to a "naive" definition of energy density 
as integration over momentum space and summation over all quantum numbers a of 
the energy (0°"^ -^°k °f created gravitons. However, the important difference is the 

appearance of the function ^^(a) which exhibits a different dependence on the scale 
factor for the zero mode compared to the KK- modes. 

Let us decompose the energy density into zero mode and KK contributions 

p = po + pKK- (43) 

Evaluating ^(a) one then obtains for the energy density of the massless zero mode 

2 f (Pk 



This is the expected behavior; the energy density of standard four-dimensional gravitons 
scales like radiation. 

In contrast, the energy density of the KK-modes at late times is found to be 

PKK = I y ^nT^l^hmnys). (45) 



which decays like 1/a^. As the universe expands, the energy density of massive gravi- 
tons on the brane is therefore rapidly diluted. The total energy density of gravitational 
waves in our universe at late times is dominated by the standard four-dimensional gravi- 
ton (massless zero mode). In the large mass limit, ninys ^ 1, n ^ 1, the KK-energy 
density can be approximated by 



Due to the factor m„ coming from the function 'S/'^, i.e. from the normalization of the 
functions ^nii.y) , in order for the summation over the KK-tower to converge, the number 
of produced gravitons ^°|^' has to decrease faster than 1/m^ for large masses and not 

just faster than 1/m^ as one might naively expect. 



4.2. Escaping of massive gravitons and localization of gravity 

As we have shown, the energy density of the KK modes scales, at late times when 
particle production has ceased, with the expansion of the universe like 

PKK-l^^ (47) 

i.e. it decays by a factor l/a^ faster than the corresponding expression for the zero mode 
graviton and behaves effectively like stiff matter. Mathematically, this difference arises 
from the distinct behavior of the functions ^%(a) and ^^,(a), n > 1, and is a direct 
consequence of the warping of the fifth dimension which affects the normalization of 
the mode functions But what is the underlying physics? As we shall discuss now, 
this scaling behavior for the KK particles has indeed a straight forward very appealing 
physical interpretation. 

First, the mass m„ is a comoving mass. The (instantaneous) 'comoving' frequency 
or energy of a KK graviton is C0n,k = + ml, with comoving wave number k. The 
physical mass of a KK mode measured by an observer on the brane with cosmic time 
dT = adt is therefore nin/a, i.e. the KK masses are redshifted with the expansion of the 
universe. This comes from the fact that m„ is the wave number corresponding to the y- 
direction with respect to the bulk time t which corresponds to conformal time r\ on the 
brane and not to physical time. It implies that the energy of KK particles on a moving 
AdS brane redshifts like that of massless particles. From this alone one would expect the 
energy density of KK-modes on the brane to decay like \/a^ (see also Appendix D of 

Now, let us define the normalized "wave function" for a graviton 

^a(^y) = ^^^. 2 rdy^>l{t,y) = 1. (48) 

y ' Jyt 

From the expansion of the gravity wave amplitude Eq. (|2TI) and the normalization 
condition it is clear that ml{t,y) gives the probability to find a graviton of mass 
for a given (fixed) time t at position y in the Z2-symmetric AdS-bulk. 

In Fig. |2] we plot the evolution of m\{t,y) under the influence of the brane motion 
with v^, = 0. 1 . For this motion, the physical brane starting at y/, ^ for t —oo moves 
towards the static brane, corresponding to a contracting universe. After a bounce, it 
moves back to the Cauchy horizon, i.e. the universe expands. The second brane is 
placed at ys = lOL and y ranges from yb{t) to y?- As it is evident from this Figure, 
is effectively localized close to the static brane, i.e. the weight of the KK-mode 
wave function lies in the region of less warping, far from the physical brane. Thus the 
probability to find a KK-mode is larger in the region with less warping. Since the effect 
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FIGURE 2. Evolution of 'i'\{t,y) = ^i{t,y)/y^ corresponding to the probability to find the first KK 
graviton at time t at the position y in the AdS-bulk. The static brane is at y^ = lOL and the maximal brane 
velocity is given by = 0.1. On the right hand panel a zoom into the bulk-region close to the moving 
brane is shown. 

of the brane motion on is hardly visible in Fig. [2l we also show the behavior of 
close to the physical brane (right hand panel). 

This shows that peaks also at the physical brane but with an amplitude roughly 
ten times smaller than the amplitude at the static brane. While the brane, coming from 
t — oo, approaches the point of closest encounter, slightly increases and peaks at 
the bounce t = where, as we shall see, the production of KK particles takes place. 
Afterwards, for ? ^ oo, when the brane is moving back towards the Cauchy horizon, the 
amplitude decreases again and so does the probability to find a KK particle at the 
position of the physical brane, i.e. in our universe. The parameter settings used in Fig. [2] 
are typical parameters which we use in the numerical simulations. However, the effect 
is illustrated much better if the second brane is closer to the moving brane. In Figured 
(left panel) we show for the same parameters as in Figure [2] but now with = L.Jn 
this case, the probability to find a KK particle on the physical brane is of the same order 
as in the region close to the second brane during times close to the bounce. However, as 
the universe expands, rapidly decreases at the position of the physical brane. 

The behavior of the KK-mode wave function suggests the following interpretation: If 
KK gravitons are created on the brane, or equivalently in our universe, they escape from 
the brane into the bulk as the brane moves back to the Cauchy horizon, i.e. when the 
universe undergoes expansion. This is the reason why the power spectrum and the energy 
density imprinted by the KK-modes on the brane decrease faster with the expansion of 
the universe than for the massless zero mode. 

The zero mode, on the other hand, is localized at the position of the moving brane. The 
profile of 00 does not depend on the extra dimension, but the zero-mode wave function 
^0 does. Its square is 




(49) 



such that on the brane (y = yb) it behaves as 

'i'lit.yb) ^ (50) 

Equation (|49l) shows that, at any time, the zero mode is localized at the position of the 
moving brane. For a better illustration we show Eq. (|49|) in Fig. [3l right panel for the 
same parameters as in the left panel. This is the "dynamical analog" of the localization 
mechanism for four-dimensional gravity discussed in IsiH]. 





FIGURE 3. Left panel: evolution of ^\{t,y) forys = L and Vf, = 0.1. Right panel: localization of four- 
dimensional gravity on a moving brane. Evolution of ^Q{t,y). Note the opposite behavior of zero mode 
and massive mode. 



This result is in contradiction with the findings of Ref. nl4\ where the authors con- 
clude that for an observer on the brane KK gravitons behave like dust with a negative 
energy density. To arrive at this result, they use Gaussian normal coordinates, 

ds^ = -N^{t,z)dt^ + Q^it,z)a^{t)5ijdx'dxj + dz^ with (51) 
Q = cosh(z/L)-7-isinh(|z|/L) = cosh(z/L) - (^y-^ - sinh(|z|/L) 

r{ty^ = ^{HLY + \ seeEq. ©. (52) 

They then argue that at low velocity, 7 — 1 , one may neglect the difference between A'^ 
and Q so that one obtains the metric 

ds^ ~ dz^ 4- e-2kl/L (_j^2 ^ ^2^^^ dijdx'dxj) . 

In this metric, the mode equation for the KK modes separates and their time evolution 
can be determined by simply solving the time part of the equation, see [14]. There is, 
however, a flaw in this argument: the above approximation is only valid sufficiently close 
to the brane (which is positioned at z = in these coordinates), but far from the brane, 
when, e.g., (7^^ — 1) sinh(|z|/L) > exp(— 2|z|/L) the above metric is no longer a good 
approximation and the difference between and Q does become important. As we have 
seen, the wave function of the KK gravitons actually is large far away from the brane 
and the time dependence enters in an important way in the normalization of the mode 
function which changes its scaling with time. 



4.3. Spectra 



In Fig. m we show the resuks of a numerical simulation for three-momentum k = 
0.01 /L, static brane position = lOL and maximal brane velocity v/, = 0. 1 . Depicted is 
the graviton number for one polarization for the zero mode and the first ten KK- 

modes as well as the evolution of the scale factor a{t) and the position of the physical 
brane yb{t)- 
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FIGURE 4. Evolution of the graviton number .yVa^k (0 for the zero mode (one polarization) and the first 
ten KK-modes for three-momentum k = 0.01 /L and Vf, = 0.1, js = IQL. 

In Fig. \5\ we show some KK spectra which we have obtained by integrating the 
equation of motion numerically. More details about the numerics and results for different 
values of the parameters can be found in Ref. LIZI . In this paper we also derive an 

10° h , ' ' ■ H 
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FIGURE 5. Final state KK graviton spectra for k = 0.001, = 100, different maximal brane velocities 
Vfo at tout = 400 for one polarization. The numerical results are compared with the analytical prediction 
(dashed line). 

analytical approximation for the spectrum which is good for KK masses m,^ < 1 . The 
numerical calculations are in very good agreement with the analytical estimates, where 
applicable. 



Integrating the zero-mode energy density over frequency with a cutoff given by the 
strong scale, k^^x — leads to the following simple result for the gravitational wave 



density parameter ||12|1 



^l/^O^yi^rad SO that Vb<0.2. (53) 

^^rad is the density parameter of the relativistic degrees of freedom at nucleosynthesis, the 
photon and three species of neutrini. The limit vt < 0.2 follows from the nucleosynthesis 
constraint which tells us that during nucleosynthesis fl^ad should not deviate by more 
than 10% from its standard value [20]. The graviton spectrum is blue with tensor spectral 
index nx = 2. Its amplitude on Hubble scales is therefore severely_suppressed and it 
leaves no detectable imprint on the cosmic microwave background [20]. 

Also the energy density of the KK modes grows like k^ for suffiently large k, 

dlogk 

and its maximum comes from the cutoff scale /cmax = We find 



a% V Prad / max \yJ \ J 

It is easy to see that low energy requires yjj < L dX all times. Therefore, to initiate a 
bounce, where yb should be close to js, we expect js^L. For typical values of the string 
scale, Ls <^L and ~ L, the above ratio is not small and back reaction of the KK 
gravitons on the geometry has to be taken into account. The ratio indicated is the one 
directly after the big bang. As time goes on the KK mode energy density dilutes faster 
than radiation and rapidly becomes subdominant. 



5. CONCLUSIONS 

In braneworld cosmology where expansion is mimicked by a brane moving through a 
warped higher dimensional spacetime, the brane motion leads to particle creation via 
the dynamical Casimir effect for all bulk modes. Here we have studied the generation of 
gravitons. 

The KK gravitons scale like stiff matter, Pkk °^ 1 /ci^, and can therefore not represent 
dark matter. In an 'ekpyrotic type' scenario with an AdSs bulk, the nucleosynthesis 
bound on gravitational waves requires Vh < 0.2. Furthermore, back reaction of KK 
gravitons on the evolution of spacetime is most probably not negligible at early times. 

In the RSII model where only one brane is present, graviton generation is negligi- 
ble tlSJ. 
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